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BaBiOa is characterized by a charge disproportionation with half of the Bi atoms possessing a 
valence 3+ and half a valence 5+. Because of selfinteraction errors, local and semi- local density 
functionals fail to describe the charge disproportionation quantitatively, yielding a too small struc- 
tural distortion and no band gap. Using hybrid functionals we obtain a satisfactory description 
of the structural, electronic, optical, and vibrational properties of BaBiC>3. The results obtained 
using GW (Green's function G and screened Coulomb potential W) based schemes on top of hybrid 
functionals, including fully selfconsistent GW calculations with vertex corrections in the dielectric 
screening, qualitatively confirm the HSE picture but a systematic overestimation of the bandgap by 
about 0.4 eV is observed. 

PACS numbers: 71.45.Lr,71.20.Nr,71.35.-y,78.20.-e,63.20.-o 



I. INTRODUCTION 

The Peierls semiconductor BaBiOa, parent compound 
of the high-T c superconductor Bai_ x K x Bi03, has long 
been of theoretical and experimental interest due to its 
distinct negative-U nature. This behavior can be under- 
stood through the concept of "forbidden valence" sug- 
gesting that in BaBiC>3 bismuth atoms appear in two 
different valences (5+ and 3+), but skip the forbidden 
formal valence in-between (4+). The "valence-skipper" 
bismuth atoms hence disproportionate, i.e. every sec- 
ond Bi atom donates all 5 valence electrons to oxy- 
gen, whereas the other Bi atom retains two 6s electrons 
and only donates the three 6p electrons to the oxygen 
atoms. Such a behavior is usually explained by the con- 
cept of a negative U which invalidates the usual positive 
Coulomb repulsion between two electrons^. The charge 
disproportionation goes in hand with significant struc- 
tural changes from the ideal cubic perovskite crystal to- 
wards a distorted monoclinic structure characterized by 
a short-ranged charge density wave (CDW) state formed 
by alternating breathing-in and breathing-out distortions 
of oxygen octahedra around incquivalent Bi 5+ (Bil) and 
Bi 3+ (Bi2) ions (Fig. [I}. As a consequence of this charge 
ordering, the formally expected metallic state for the 
cubic perovskite BaBi 4+ 03 is replaced by an insulating 
regime characterized by a large direct CDW optical re- 
sponse Ed « 2.0 eV— and an indirect optical transition 
E;. Although BaBi03 has been the subject of numerous 
optical investigations, dissenting opinions were reported 
on the value of the indirect gap Ei, for which conflicting 
measurements are reported in literature, ranging from 
0.2 e\£ and 1.1 eV4 

Despite the apparent absence of strong electron local- 
ization usually associated with d and / electrons, stan- 
dard semilocal density functionals exhibit severe prob- 
lems in describing the basic properties of this mate- 
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FIG. 1: (Color online) Self-consistent-GW charge density 
plots showing (a) the breathing distortion in the (100) plane 
and (b) the tilting instability projected in the (001) plane. 
Bi 5+ (Bil) and Bi 3+ (Bi2) indicate the two inequivalent bis- 
muth ions, whereas the arrows indicate the breathing dis- 
placement of the oxygen atoms, (c) Schematic view of the 
BaBiOa unit cell showing the BiOe tilted octahedra. Color 
coding: darker (red) areas indicates a high value of the charge 
density. 



rial, yielding much too small or even negative band 
gaps^£, thus preventing a detailed description of the elec- 
tronic properties. In a recent paper, we have explained 
the polaron-mediated insulator-to-metal transition ex- 
perimentally observed in K-doped BaBiC>3 by means of 
hybrid density functional theory. We have shown that 
semilocal functionals are incapable (i) to disclose the 
underlying physics of this type of electronic transition 
and, more generally, (ii) to account for "long range" 
Peierls-like distortions in doped BaBiO^. The objective 
of the present study is to apply the screened exchange 
hybrid density functional theory (DFT) Heyd-Scuseria- 
Ernzerhof (HSE) functional 7 - to the ground state struc- 
tural, electronic and vibrational properties of BaBiC>3, 
and to investigate to which extent many-body effects may 
alter the ground state and quasi-particle (QP) properties 
of BaBiC>3. We address the latter issue by adopting the 
GW approach. These calculations essentially confirm the 
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picture emerging from the hybrid functional. Details of 
our computational approach are described in the method- 
ological section, whereas the results are presented and 
discussed in Sec lIIIl 



II. METHODOLOGY 

We have used the projector augmented wave method^ 
based Vienna ab initio simulation package (VASP)^ em- 
ploying the generalized gradient approximation scheme 
of Perdew et alM (PBE), the HSE hybrid functional^ 
and the GW formalismi2.Il The HSE functional is 
constructed by mixing 25% of the exact Hartee-Fock 
(HF) exchange with 75% of the Perdew-Burke-Ernzcrhof 
(PBE) exchange functional 7 . The resulting expression 
for the exchange (x) and correlation (c) energy is 



^xc ^x 4 33 



HF,sr,/i 



+ £ C PBE (1) 



in which the parameter fi controls the decomposition of 
the Coulomb kernel into short-range (sr) and long-range 
contributions to the exchange. E'HF.sr,^ j s the exact ex- 
change screened at very long range, and ££ BE > sr >' J is 
the corresponding short range exchange functional con- 
structed using the PBE exchange hole. According to the 
HSE06 recipe, we set fj, = 0.2 A" 1 ^ This implies that 
non-local exchange is rather long ranged extending over 
several nearest neighbor shells. The HSE approach has 
been used to calculate the starting orbitals for different 
GW-based schemes of increasing complexity and preci- 
sion in the evaluation of the self-energy E = iGW. The 
screened Coulomb kernel W — e _1 w which enters in the 
previous expression requires the computation of the fre- 
quency dependent dielectric matrix e _1 = 1 + vx, where 
the polarizability x is given by the Dyson equation 



X = [1 - Xo(w + fxc)} 1 Xo 



(2) 



Here, xo indicates the independent particle polarizability, 
whereas f xc represents the effective non-local exchange 
correlation kernel which describes the many-body inter- 
action between electrons and holes (vertex corrections). 
We have adopted three different GW strategies: 

(i) The widely used single shot (i.e. non-self- 
consistent) Go Wo approximation, where Go and Wo 
are calculated using HSE eigenvalues and eigenfunctions 
within the so called random phase approximation (RPA) 
(i.e. excluding any vertex corrections, /z C =0)J£ The 
VASP implementation is detailed in Refs. [TH and [Til- 
Using HSE orbitals, this approach usually slightly over- 
estimates the band gap, but is very efficient ji 7 - 

(ii) A partially self-consistent procedure consisting 
of an update of the eigenvalues in the Green's func- 
tion G combined with HSE screening properties Wo 
(GWo). Within this scheme we have included electro- 
static electron-hole interactions in the computation of 
W using time-dependent HSE (TD-HSE) through an ef- 
fective non-local frequency- independent kernel f xc . This 



specific GW scheme, usually denoted as GWj 01 " (test- 
charge/test-charge GWo), is believed to remedy most of 
the problems of single shot Go Wo and is considered to 
be a convenient (less time-consuming) and trustable (ex- 
perimental precision) alternative to fully self-consistent 
GW. The approach was found to yield excellent QP 
band gaps in simple sp bonded semiconductors, if and 
only if the electrostatic interactions between electrons 
and holes were taken into account ^ The electrostatic 
electron-hole interaction is included using an effective 
many-electron exchange correlation kernel constructed 
using the HSE functional (Nano-quanta kernel deter- 
mined to mimic TD-HSE, for details we refer to Refs. 

(iii) Finally, we employ the self-consistent quasiparti- 
cle GW (scQPGW TCTC ), following the original ideas of 
Faleev, Schilgaarde, and Kotan i 23 i 24 , with the inclusion 
of vertex corrections adopting the formalism of Reining 
et a/i 19 ' 20 i 21 ' 22 [Nano-quanta kernel constructed to mimic 
the Bethe-Salpeter equation (BSE)]. This approach is 
rather time-consuming but entirely free of any depen- 
dence on the starting orbitals. For sp bonded materials, 
it yields band gaps within a few percent of experiment, 
if the electrostatic electron-hole interaction is included 
using a many-electron exchange correlation kernel f X c^ 

The simulated unit cell, containing 10 atoms, has been 
modeled using a 4x4x4 fc-point grid and an energy cutoff 
of 400 eV, reduced to 300 eV within GW calculations. 
To calculate the frequency-dependent dielectric matrix 
the inclusion of 400 empty bands was sufficient to obtain 
well converged screening properties and QP eigenvalues. 



III. RESULTS AND DISCUSSION 

A. Structure 

The low temperature phase of BaBiOs is monoclinic 
and derives from the primitive cubic perovskite high tem- 
perature cell by simultaneous breathing (6) and tilting 
(4>) distortions of the octahedra 26 (Fig. [T|). Unlike pre- 
vious studies, we have performed a full structural relax- 
ation by allowing for both b and 4> distortions as well 
as for volume (V) and shape (/3, the monoclinic angle) 
optimization. The resulting minimized geometries listed 
in Tab. U reveal that the structural properties are sensi- 
tive to the applied functional and are generally better re- 
produced within HSE. In particular, the PBE functional 
seriously overestimates the volume (« 4 %), which is in 
fact well described using HSE. Also, the breathing dis- 
placement, which is responsible for the opening of the 
CDW gapl is better described using the hybrid func- 
tional. The minimum HSE Bi-0 bond length splitting b 
= 0.09 A matches exactly the experimental value. The 
volume and, to a lesser extent, the breathing displace- 
ment does not improve by using a revised version of the 
PBE functional designed to yield improved equilibrium 
volumes, the so called PBEsol functional 2 ! (see Tab. [I]). 
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TABLE I: Collection of calculated quantities along with the available experimental data. The abbreviations V (volume), b 
(breathing distortion), (j> (tilting distortion), j3 (monoclinic angle), Ed (direct gap), Ei (indirect gap) andtoo (macroscopic 
dielectric constant) are described in the text and, partially, in Fig. [T](c). a RefEE 6 RefH c Ref[3 d RefJH e RefSl / RefQ], 9 



Refl30j, h Ref|31|, 





PBE 


PBEsol 


HSE 


G Wo GWj CTC 
(This work) 


scQPGW TCTC 


LMTO-LDA 
(Ref.4) 


PAW-LDA 
(Ref.5) 


Expt. 


v (A 3 ) 

6(A) 


85.76 
0.07 
12.1 

90.16 


85.94 
0.075 
12.2 
90.14 


82.10 
0.09 
11.9 

90.24 






80.41 
0.11 
9.6 


81.04 
0.04 


82.21 a 
0.085 a 
11.2 a 
90.17 a 


E d (eV) 
Ei (eV) 


1.22 
0.0 


1.32 
0.0 


2.07 
0.84 


2.62 2.39 
1.32 1.15 


2.45 
1.28 


1.1 
0.1 


0.0 


2.0 b 
0.2 c , 0.5 d , 0.8 


Coo 


12.2 




6.0 


7.0 


7.2 






5 9 , 7 h 



B. Electronic properties 

The approximations made in the exchange-correlation 
term, have an even more drastic impact on the predicted 
electronic and dielectric properties. In line with past find- 
ings, PBE predicts a zero-gap system (E BBE = 0.0 eV) 
with a much to small direct band gap Ed (E PBE = 1.22 
eV) . Although the origin of this failure has often been re- 
lated to the deficient description of the structural insta- 
bilities at the PBE level, we found that the values E PBE 
and E BBE computed adopting the HSE and experimen- 
tal geometries are only 0.05 meV larger than those cal- 
culated at the PBE structural minimum, and therefore 
still in conflict with the experimental measurements. 

In order to capture the insulating character of BaBiC>3 
we have adopted the HSE functional and the GW method 
at the three levels described in the methodological sec- 
tion. The computed band structure is shown in FigJ^l 
where we plot the most relevant bands near the Fermi 
level (E F =0) for PBE, HSE, G W n , GWj CTC and 
scQPGW TCTC . We observe that the pragmatic inclusion 
of 25% non-local Hartree-Fock exchange within the HSE 
approach shifts apart the bands near Ep and increases 
both Ei (0.84 eV) and E d (2.07 eV). These values are in 
good agreement with the available measurements. 

The GoWo approximation evaluated using HSE or- 
bitals further increases the band gap between the oc- 
cupied and empty states resulting in significantly larger 
quasiparticle gaps than HSE (see Tab. HJ). The update 
of the eigenvalues in G and, in particular, the inclusion 
of e-h (electron-hole) interactions (GWj CTC ) leads to a 
slight reduction of the quasiparticle band gaps (« 4%). 
We note that the former approach (Go Wo) is known 
to overestimate the band gaps for s-p bonded semicon- 
ductors and insulators using HSE orbitals by typically 
10 whereas the latter approach (GWj 01,0 ) yields 



very good values almost on par with fully self-consistent 
scQPGW TCTC calculation o 18 ! 25 . Indeed, in this specific 
material, the scQPGW TCTC band structure is very close 
to the GWj band structure, although the indirect 
band gap increases slightly from 1.15 eV to 1.28 eV. Un- 
fortunately, the experimental uncertainty in Ei does not 
permit to ascertain which method yields the most accu- 
rate results, but the GW values are generally somewhat 
larger than the maximum values determined experimen- 
tally. Since the direct gap Ed has been measured using 
optical techniques, it is possibly impaired by excitonic 
effects (see below), but overall it seems that GW does 
yield a slightly too large band gaps for BaBiOa. However, 
compared to photoemission data—, the calculated total 
density of states (shown on the right side of Fig[2]) show 
that GWj 01,0 and scQPGW TCTC yield a main valence 
band peak (w -3.4 eV) in much better agreement with 
experiment (-3.2 eV) than PBE (-1.9 eV) and HSE (-2.8 
eV). Once more, we underline the need of more detailed 
experimental data to unequivocally probe the predictive 
power of GW-type and HSE calculations. It should also 
be stressed that GW-like approaches, though employed 
for more than 20 years, have been applied so far mainly to 
prototypical semiconductors and insulators (GaS, ZnO, 
CU2O etc..) but, to our knowledge, their performance 
on more complex systems such as BaBiOa, has not yet 
been assessed in literature. A comparison with accurate 
experimental measurements is therefore essential. 

Let us proceed to shed light on the nature of the CDW 
gap Ei. It is accepted that Ei is formed by the splitting 
of the conduction Bi(s)-0(p) anti-bonding orbitals into 
two subbands, but the electronic charge redistribution 
has not been discussed in details. Figure [3] (a) reports 
the scQPGW TCTC Z-projected density of states (DOS) 
on Bi and O sites which should be interpreted in con- 
junction with the band-projected charge (Tab. [TT|. The 
formation of the CDW lowers the symmetry of the crys- 
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FIG. 2: (Color online) Calculated (left) band structure of monoclinic BaBi0 3 within PBE, HSE G W , GWq and 
scQPGW and (right) corresponding total density of states. Ei is opened between the A and R points of our modeled 

supercell, which correspond to the W and L points of the simple monoclinic cell. P refer to the measured main photoemission 
peak of the valence band^. 
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FIG. 3: (Color online) Calculated scQPGW TCTC (a) partial 
DOS decomposed over Bil, Bi2 and O (apical and planar 
Oxygen yield identical DOS) and (b) DOS difference between 
Bil and Bi2 decomposed into / quantum numbers. 



TABLE II: scQPGW TCTC , HSE and PBE band-projected 
charge (in partial number of electrons e~) on Bil and Bi2. 
Charges have been determined inside the PAW spheres. 



P 



Bil 
Bi2 

Bil-Bi2 



Bil 
Bi2 

Bil-Bi2 



1.280 
1.642 
-0.36 



1.318 
1.694 
-0.38 



scQPGW 
1.199 
0.871 
+0.33 

HSE 
1.176 
0.830 
+0.35 



TCTC 



10.254 
10.065 
+0.19 



10.252 
10.063 
+0.19 



tal thus inducing a global rearrangement of the charge 
around the two inequivalent Bi sites (see Fig. [T]) . We find 
that Ei is opened between occupied Bi2(s) states (va- 
lence band maximum) and empty THl(s)-Bi2(p x p y ) or- 
bitals (conduction band minimum). This orbital-ordered 
insulating state is indeed more complex than the gener- 
ally accepted picture invoking a simple Bi(s)-0(p) gap 
and arises from the elaborate nature of the charge dis- 
proportionation, resulting in the electronic states shown 
in Figured (b). 

Previous calculations^^ argued that the static valence 
charge difference 5 between Bil and Bi2 is not significant. 



Bil 
Bi2 

Bil-Bi2 



1.371 
1.623 
-0.25 



PBE 
1.085 
0.863 
+0.22 



10.181 
10.056 
+0.13 



The calculated S is indeed rather small (5pbe 
$hse = 0.16e~, ScQPGW 1 



0.16e") but 



= O.lOe", 
it is the 



result of a substantial transfer involving orbitals with 
different angular moments. HSE and GW provide es- 
sentially the same picture, whereas PBE underestimates 
significantly the ^-decomposed charge transfer by about 
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30% (see Table |TTJ) . In more detail, the flow of charge 
induced by the charge disproportionation between Bil 
and Bi2 involves a s-like charge transfer from Bil to Bi2 
(0.36e~) which is compensated by a backtransfer of p 
(0.33e~) and d (0.16e~) charge mostly localized in the 
energy range from -8 and -4 eV. This backtransfer can 
be understood by realizing that the O atoms are much 
closer to Bil. As a result, the tails of the 0(p) orbitals 
overlap with the Bil spheres and are picked up as states 
with higher angular momentum specifically p and d like 
character. 



C. Optical and dielectric properties 




son 



FIG. 5: (Color online) Comparison between calculated and 
measured^ infrared reflectivity. 
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FIG. 4: (Color online) The imaginary parts of the di- 
electric function (e 2 ) are shown for RPA-HSE, TD-HSE 
and scQPGW TCTC along with the experimental curve from 
Ref[H. 

In addition to the formation of a band gap, the 
CDW also causes a strong optical resonance at w 2 
eV. In Fig|4] we compare the experimental imaginary 
part of the dielectric function with the theoretical 
spectrum computed within different approximations: (i) 
the random phase approximation (RPA) using HSE or- 
bitals, (ii) TD-HSE (as described in Ref. 0j) and (iii) 
TD-scQPGW TCTC . The first approach (i) neglects any 
electrostatic interactions between electrons and holes. 
The second approach (ii) is similar to TD Hartree-Fock 
but the electrostatic interaction between electrons and 
holes is only 1/4 — more precisely it is described by the 
sr-exchange kernel in HSE — whereas in the final ap- 
proach (iii) the interaction between electrons and holes 
is described fully ab-initio using the W determined in 
the scQPGW calculations. For (i) and (ii) HSE or- 
bitals and one-electron energies are used, whereas (iii) 
is based on the scQPGW one electron energies and or- 
bitals, and the results of (iii) are usually equivalent to 
Bethe-Salpeter calculations without the Tamm-Dankof 
approximation j 18 i 20 



Considering the complexity of the system the agree- 
ment is remarkably good with regard to the position and 
the width of the two dominant peaks namely the CDW 
peak at 2 eV and the excitation at 8.0 eV, which can 
be described as the transition from the main peak of the 
valence band (-3.4 eV) to the broad group of states start- 
ing at about 4 eV above Ep. The overestimation of the 
intensity of the CDW peak is most likely related to a 
too coarse k-point sampling. Unfortunately, in this fairly 
complex system, calculations with a denser k-point grid 
are very time consuming, and we do not expect that they 
will change the results qualitatively. The first important 
observation is that there is a difference between the RPA- 
HSE and TD-HSE calculations. Excitonic effects (e-h in- 
teraction) shift the CDW peak down by approximately 
0.3-0.4 eV. The scQPGW calculations, which also include 
excitonic effects, predict a blue shifted CDW excitation ss 
0.2-0.4 eV above the experimental peak. Unfortunately, 
this confirms that the scQPGW overestimates the direct 
band gap somewhat, as was already hinted at in the sec- 
tion on the electronic properties. 

Finally, we calculated the ion-clamped (high fre- 
quency) macroscopic dielectric constant (Tab]!]). The 
HSE too was computed adopting the perturbation expan- 
sion after discretization (PEAD) method 35 . Within the 
PEAD method we applied a finite field in the supercell 
and extracted the static dielectric constant from the re- 
sultant macroscopic field. This approach includes local 
field effects (electron- hole interactions). The so obtained 
macroscopic dielectric constant (e 00 ffS£ '=6) is in good 
agreement with the experimental estimations (5 and 7, 
see TabJU extracted through a fit of the reflectivity spec- 
tra. The GW and scQPGW values, 7.0 and 7.2 respec- 
tively, though still in good agreement with experiment 
are about 15% larger than the HSE values. As expected 
from the serious underestimation of the gap, TD-PBE 
fails dramatically in predicting eoo (e oc PBE =12.2). 
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D. Phonons and infrared reflectivity 

We conclude our analysis by exploring the effect of non- 
local exchange on the zone-centered phonon frequencies 
and on the oscillator strengths of the infrared reflectiv- 
ity (IR) active modes. Since within the GW method, 
the calculation of phonon frequencies is presently out of 
reach, we focus on PBE and HSE only. The evaluation of 
the Born effective charge (BEC) Z* a/3 , required for the 
calculation of the IR spectra, were performed adopting 
the PEAD (HSE) and linear response (PBE) methods. 
Excluding the three acoustic modes (which are zero at 
r), the monoclinic phase (containing 10 atoms per unit 
cell) possesses 27 optical modes with the symmetry 

7A g + hB g + 6A U + 9B U 

. The gerade modes are Raman active and the ungerade 
modes are infrared active. In Fig. Owe compare our cal- 
culated HSE IR spectra with the experimental reflectivity 
data2£. The results are remarkable: the inclusion of non- 
local exchange is essential to account for a correct predic- 
tion of the phonon modes with large oscillator strength. 
The three highest modes at 157, 256 and 463 cm -1 re- 
produce very well the corresponding experimental values, 
whereas the lowest mode at 119 cm" 1 is shifted upward 
with respect to the measured one by w 20 cm -1 . In ad- 
dition, we find that the highest frequency phonon (Ra- 
man active) corresponding to the BiOe breathing mode 
responsible for the CDW instability, experiences a huge 
shift from PBE (453 cm- 1 ) to HSE (562 cm" 1 ). The 
HSE value coincides with the value provided by Raman 
scattering analysis^, 570 cm -1 . The problems of PBE 
can be attributed primarily to volume effects: Due to the 
large PBE volume the PBE phonon vibrations are located 
at lower frequencies. In addition to this, we note that 
the PBE inadequate description of the oxygen breath- 
ing distortions is also reflected in the phonon spectra: 
The frequency downshift is larger (45-60 cm -1 ) for the 
oxygen-based modes lying at higher frequency than for 
the lower modes (20-30 cm -1 ) emerging from the softer 
Bi and Ba vibrations which are actually very well de- 
scribe at the PBE level. A similar behaviour has been 
observed in ZnO^ 

Finally, we list in Tab. [ITTI the computed BEC Z* aj} . 
We find that the BEC tensor is diagonal and reduces to a 
scalar matrix for Ba and Bi (i.e. Z* Ba/mil = Z* Ba/Bi22 = 

^Ba/Bi 33) an( i t° tw° components Zt and Z± for O with 
respect to the Bi-0 bond. In agreement with prece- 
dent estimations of Thonhauser and Rabe^. we observe 
a significant BEC disproportionation between Bil and 
B2 (Z Bi2 — Z Bil ), substantially larger than the corre- 
sponding static charge transfer 5 discussed above (see 
Tab. mj). The incorporation of non-local exchange ef- 
fects through HSE reduces the dynamical charge dispro- 
portionation from 1.66 |e"| (PBE) to 0.97 |e"| (HSE) 
and leaves the Ba and O BEC's unchanged. The aver- 
age BEC of oxygen -3.13 |e~| equals the experimentally 
derived value of -3.2 |e - |2i. 



TABLE III: Calculated HSE and PBE Born effective charges 
Z* a/3 (with i atom index) in unit of \e~\. and Zj_ refers to 
the components parallel and perpendicular to the Bi-O bond. 
The LDA values are taken from Ref. 0. 



Bil 



5.01 



Bi2 



Ba 



Z\ 



Zi }C 

z\ 



z\ 



4.65 



4.78 



HSE 
5.98 2.71 



PBE 
6.31 2.71 



LDA r ' 
6.22 2.75 



O 



-4.39 
-1.88 



-4.36 
-1.89 



-4.55 
-1.85 



IV. SUMMARY 

In summary, by means of beyond-DFT methods we 
have computed ground and excited state properties of 
BaBi03. On the basis of hybrid functionals and fully ab 
initio GW methods, we have shown that the inclusion of 
non-local exchange effects is essential to reproduce and 
understand the complex physical properties of BaBiC^. 
Our calculated structural distortions (oxygen instabili- 
ties and cubic-to-monoclinic transition), optical excita- 
tions (direct /indirect band gap and optical spectrum), 
dielectric (dielectric constant) and vibrational properties 
(phonon frequencies, infrared reflectivity and Born ef- 
fective charges) reproduce well the available experimen- 
tal findings. Overall, our study suggests that the GW 
approaches tend to overestimate the band gap, which 
seems to be better predicted within the screened hy- 
brid HSE functional. Although the inclusion of an ef- 
fective non-local exchange-correlation kernel f xc at TD- 
HSE (GWj CTC ) and BSE-GW (scQPGW TCTC ) level im- 
proves the single-shot GqWq description, the resulting 
band gap is still « 15% (E rf ) and 40% (E<) larger than 
the corresponding bare HSE value. The source of this 
error as far as it can not be attributed to uncertainties 
in the experimental results likely resides in the neglect of 
many-body vertex corrections in the self-energy £, which 
are entirely neglected in this work, or any other GW 
implementation. Though scQPGW TCTC has proven to 
be excellent for the prediction of the band gaps of typi- 
cal semiconductor, insulator and noble gas solids it was 
never applied to more complex oxides such as BaBiOa. 
Given the lack of accurate experimental information on 
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the electronic properties of BaBiC>3 (especially for the in- 
direct gap Ej) the assessment of these GW techniques to 
less critical oxides for which well established experimen- 
tal data exists is highly required. 

Concluding, our research on BaBiC"3 clarifies the se- 
rious drawbacks of standard DFT and provide a key to 
understand and design similar charge-ordered materials. 



We hope that our study can boost for new experiments 
as well as for thorough theoretical and computational ef- 
forts aiming to straighten the issues still under debate. 
Furthermore, we have shown that the hybrid functional 
HSE yields a description which is at least on par with 
sophisticated man-body techniques. 
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